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Abstract 

We study the decays B K^*^£'^i~' and Bg £^i~ {£ = e,/i) in the minimal 
supersymmetric standard model with R-parity violation (RPV). From the recent mea- 
^ • surements of their branching ratios, we have derived new upper bounds on the rele- 

vant RPV coupling products, which are stronger than the existing ones. Using the con- 
^ ■ strained parameter space, we predict the RPV effects on the forward-backward asymme- 

. 5^ ; tries Afb{B K^*k+£-) and the branching ratios B{Bs i+£~). Our results of the 

forward-backward asymmetries agree with the recent experiment data. It is also found 
that B{Bs — > £'^£~) could be enhanced several orders by the RPV sneutrino exchange. 
The RPV effects on the dilepton invariant mass spectra of B ^ K^*k+£- and the nor- 
malized Afb{B ^ K(-*k+£~) are studied in detail. Our results could be used to probe 
RPV effects and will correlate with searches for direct RPV signals at LHC. 
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1 Introduction 



Flavor changing neutral current (FCNC) b ^ s processes are forbidden at the tree level in 
the standard model (SM), which proceed at a low rate via penguin or box diagrams. If addi- 
tional diagrams with non-SM particles contribute, their rates as well as other properties will be 
modified. This feature make FCNC processes powerful means to probe new physics indirectly. 
The recent experimental measurements of 5 ^ K^*^i'^i~ decays [H [21 El H] agree with the SM 
predictions within their error bars, therefore, these measurements will afford an opportunity to 
constrain new physics scenarios beyond the SM. 

Semileptonic rare decays B K^*^i~^i~ have been extensively studied previously. The 
dominant perturbative SM contribution had been evaluated years ago[S], and later QCD cor- 
rections have been provided P,[7l[H]- Cf^l/ml) corrections have been first calculated in Ref.p] 
and then in Refs.[TOl ITT]. Long distance contributions can have different origins according to 
the value of the dilepton invariant mass. The contributions of charmonium resonances to these 
decays by means of Vector Meson Dominance (VMD) have been studied carefully [TUl \n\ [T ^ [T^ . 
Far from the resonance region, instead, cc long-distance effects are investigated using a heavy 
quark expansion in inverse powers of the charm-quark mass ((9(l/m^) corrections) [15]. Anal- 
yses of new physics contribution have been performed in different models, for example, the 
two-Higgs doublet model[16] , the supersymmetric (SUSY) models [ITIIIH], the SUSY SO(IO) 
grand unification theory [19] and the top quark two-Higgs doublet model [20] . 

The effects of RPV SUSY in B meson decays have been extensively investigated in the 
literature [211 [22] • The decays B K^*H~^£~ and Bg are all induced at the parton 

level by 6 ^ si~^i~ process, and they involve the same set of the RPV coupling products. In 
this paper we will study the decays B K^*U^i- and 5, i+i' in the RPV SUSY model. 
Using the recent experimental data, we will obtain the new upper limits on the relevant RPV 
coupling products. Then we will use the constrained regions of the parameters to examine 
the RPV effects on the branching ratios of Bs decays and the forward-backward 

asymmetries (Afb) oi B K^*H~^£~ . In addition, we will compare the SM predictions with the 
RPV predictions about dilepton invariant mass spectra and the normalized forward-backward 
asymmetries in B ^ K^*^i~^i~ decays. 

The paper is arranged as follows. In Section 2, we introduce the effective Hamiltonian and 
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calculate the expressions for B K^*^£~^£~ and £^£~ processes in the RPV SUSY. In 

Section 3, we tabulate the theoretical inputs and deal with the numerical results. We display 
the constrained parameter spaces which satisfy all the available experimental data, and then 
we use the constrained parameter spaces to predict the RPV effects on Afb{B -> K^*H+£-) 
and B{Bs £~^£~), which have not been well measured yet. We also show the RPV effects 
on dilepton invariant mass spectra and the normalized forward-backward asymmetries in i? — >• 
}{(*)£+£- decays. Section 4 contains our summary and conclusion. 



2 The theoretical frame for B K^*k+i- and Bs 

2.1 The decay branching ratios in the SM 
2.1.1 The semileptonic decays B K^*^£+£- 

In the SM, at the quark level, the rare semileptonic decays h s£^£~ can be described by the 
effective Hamiltonian 

^ 10 

Hf/f (6 - s£n-) = -^v*v,, c.if^n, (1) 

where the operator base Oi is given in [7j. We will use the Wilson coefficients Cj(yu) calculated 
in the naive dimensional regularization (NDR) scheme [7j, and the long-distance resonance 
effects on C^^^ {^) given in Refs.|12l [T^ . The Hamiltonian leads to the following free quark 
decay amplitude: 

M^'\h ^ s£+£-) = 9l^v:yt,{clf\s^,PLh){£Y£) + C,o{s^,PMh'ld) 

-2m,C'/^ (^sta^XPnb^ {£Yi)}, (2) 

with Pl,r = (1 =F 75)/2,s = and q = p+ + p„ {p± the four-momenta of the leptons). In 
our following calculations, we take ms/rrif, = 0, but keep the lepton masses. The hat denotes 
normalization in terms of the B-meson mass, m^, e.g. s = s/m\, rhq = mg/niB- 

Exclusive decays B K^*^£~^£~ are described in terms of matrix elements of the quark 
operators in Eq.([2]) over meson states, which can be parameterized by the form factors. It is 
worth noting the form factors involving the B — > K^*^ transitions have been updated recently 
in EHl. 
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Using Eq.([2]), one can get the amplitudes of exclusive B K^*H~^i decays 



where for B K£^i~, 



TT 



(3) 

(4) 
(5) 

Ti^ = A{s)e^p^pe*Pp%pfl,, -iB{s)el + iC{s){e* ■pb)p,+ (6) 
= E{s)e,,^pe*''p'ip^j,, - tF{s)e; + ^G(5)(e* ■ pb)p, + ^i/(5)(e* ■ pB)q„ (7) 

with p = pb +PkM- Note that, using the equation of motion for lepton fields, the terms in 
vanish, and those in become suppressed by one power of the lepton mass. 



and for B K*i+t 



The 



auxiliary functions in 7^^ and 7^^ are defined as 



i + m/< 



B'(s 



C'{s) 



s 

= Cio/+(S), 



D'{s) = Cio/-(5), 

A{s) = —^C',"{s)V{s) + ^C^/fT,{s), 

s 



B{s) = {l+niK' 
C{s) = ^ 



{l-mK*)Cl"{s)A2{l 



m = 



1 — mfj-* [ \ 

C'^"{s)[{l + mK*)A,{s) - (1 - mi,.)^2(s) - 2mK'Ao{s 



1 - mx* 



1 



1 + rriK* 
F{s) = {l+mK*)CioA^is 

Gis) = \ CioA2(g), 
1 + mi^-. 



//(s) = lCio\{l + mK*)Ai{s)-{l-mK*)A2{s)~2mK*Ao{s] 

S L 



(8) 

(9) 
(10) 

(11) 
(12) 

(13) 
(14) 

(15) 
(16) 
(17) 
(18) 

(19) 



It's noted that the inclusion of the full s-quark mass dependence in the above formulae can 
be done by substituting nib ^ nib + fns into all terms proportional to Cj^^Ti and Cj^-^ fx and 
mb-^nib- in C^/%s, since ~ sa.^inib + m.) + (m, - m,H]q'^b. 
The kinematic variables (s,u) are chosen to be 



u 



f = {P++P-) 



{PB-P~f - {pB-p+f: 



(20) 
(21) 



which are bounded as 



(2m^)2< s <(l-m^w)' 



-uis] < u < uis] 



(22) 
(23) 



with rhi = mi/mB and 



u{s) = ^A(l-4^) 
A = A(l,m^(,),s) 



1 + m^(.) -2s- 2m^(.)(l + s) 



(24) 



(25) 



The variable u corresponds to 6, the angle between the momentum of the B-meson and the lep- 
ton in the dilepton center-of-mass system (CMS) frame, through the relation u = —u{s)cos6 
[14] . Keeping the lepton mass, we find the double differential decay branching ratios and 
B^* for the decays B — Ki~^i~ and B K*i~^i~ , respectively, as 



dsdu 



Gl^alnv" 



x|(|Af + |Cf)(A-n2) 



dsdu 



+ \C'\'Ami{2 + 2mi- - s) + Re{C'D'*)87fii(l - mj^) + \D'YAnijs 



(26) 



Gl^alnv" 



tb- 



X 



UP 



E\ 



(s(A + u^) + Ani]\) + ^(s(A + u^) - Arhlx) 



1 



4m|^, 



fin A - + 8m^*(s + 2m^2) ) + \F\^[\ - + 8m^*(s - 4m 



-25^ 



Re{BE*) + i?e(AF*) 



A 



1 



Re{BC*){l - ml, - s){X - u^) 
- Re{FG*) ((1 - m|.* - s){X - u^) + 



rrie 



-A 



K* 



Re{FH*) - Re{GH*){l - m^.] 



+ \H 



2 



■sA . 



m 



(27) 



K* 



Our results of the double differential decay branching ratios are consistent with the ones in 
Reim. 



2.1.2 The pure leptonic decays Bg i^i 

In the SM, the effective Hamiltonian have been given by 



VlVtbY{xt){-sh)v-Amv-A + h.c. 



(28) 



where Xt 



,mt = mt{mt), {sb)v-A = ■57/^(1 — Ts)^. The pure leptonic decay amplitudes 



can be written as 



M^'^iB, _ i+r) = hsM {i ^b{i- 



(29) 



with 



h 



Up 



SM 



v/2 27rsin^6'iy 
Then we can get the branching ratios for Bg —>■ 



(30) 



1 - Amj 



71 V47rsin^^ 



w 



M''\B, 

2 



flminiBjl - 4mj \V;M'y\^t). (31) 



2.2 The decay amplitudes in RPV SUSY 

In the most general superpotential of the minimal supersymmetric Standard Model (MSSM), 
the RPV superpotential is given by [25] 



(32) 
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where L and Q are the SU(2) doublet lepton and quark superfields, respectively. E'^, If^ and 
D'^ are the singlet superfields, while i, j and k are generation indices and c denotes a charge 
conjugate field. 

The bilinear RPV superpotential terms fiiLiH^ can be rotated away by suitable redefining 
the lepton and Higgs superfields [26]. However, the rotation will generate a soft SUSY breaking 
bilinear term which would affect our calculation through loop level. However, the processes 
discussed in this paper could be induced by tree-level RPV couplings, so that we would neglect 
sub-leading RPV loop contributions in this study. 

The A and A' couplings in Eq. ( l32l) break the lepton number, while the A" couplings break the 
baryon number. There are 27 A^^^ couplings, 9 A^fc and 9 Af^-^ couplings. Xiij]k are antisymmetric 
with respect to their first two indices, and A^'j^-^] are antisymmetric with j and k. 




K,K* e 




K,K* 



(3) (b) 

Figure 1: The RPV contributions to i? ^ K^*H^i~ due to sneutrino and squark exchange. 





(a) 



Figure 2: The RPV contributions to Bg — > i^i due to sneutrino and squark exchange. 




From Eq. fl32l) . we can obtain the relevant four fermion effective Hamiltonian for b — > sii 
process due to the squarks and sneutrinos exchanges 



1 V- 



A •fcA'* 

U ^ imn (^d^p^d^^ [ikPl^j] 



m 



(33) 



The RPV feynman diagrams for B K^*^i^i~ and Bs are displayed in Fig{T] and 

FiglSl respectively. 

From Eq.( l33|) . we can obtain the RPV decay amplitude for B Ki'^i' 

M^^ [B ^ Ktn = As(4(^ij+ ^k){1 - 75)^,) + Ap(4(l - 75)^,) + A',(4(l + 75)^,),(34) 
with 



■-f 



\* \/ 2 2 



(35) 
(36) 
(37) 



nib - ms 

and nib and nis the quark's running masses at the scale rUb. 
For B ^ K*e+e~, the RPV amplitude is 

M"^^ {B ^ K*i+r) = - 75)^i) + ^^^(4(1 - 75)^j) + ^^',(4(1 + 75)^,), (3^ 

where 



E 



^ 8m? 



B-*K* , 
' 



S , 1 



2 + nis 

K* 



B 



2 nib + 



"5J xi 2 

X'inig 



and the auxiliary functions are given as 



m 

J{s) 
Kis) 

m 



E 

i 

E 

i 

E 

i 

E 



^j^t3^fc^2 
^j^t3^fc^2 
^ji3^fct2 
^ji3"^fct2 



2r 



-m. 



-(mB + m;^O^f"'''*(S) 



A 



B^K' 



-m 



B 



rriB + rriK* 
'2.niK* ( ,B-^K 



■{Ai 



A 



B~*K* 



(39) 
(40) 

(41) 



(42) 
(43) 
(44) 
(45) 



For Bs i+£-, the RPV amplitude is 
M"^^ {Bs ^ tr) = Au{ik ^b{1 - 75)^,) + Ap(4(l - 75)^,) + A',(4(l + 75)^,), 



(46) 



with 



As 

a:. 



E 



-if'Bsf^Bs), 



and 



-Is 



The RPV couphngs can be complex in general, we write their products as 



where the RPV coupling constant A G {A, A'}, and <^^^ is the RPV weak phase. 



(47) 
(48) 
(49) 



(50) 



2.3 The branching ratios with RPV contributions 

With these formulae in Sec. 2. 1-2.2, we can obtain the total double differential decay branching 
ratios of the decays B K^*H~^i~, 



(51) 

dsdu dsdu ' dsdu ' dsdu ' dsdii 
Since we will only consider one RPV coupling product contributes at one time, we have neglected 
the interferences between different RPV coupling products, but kept their interferences with 
the SM amplitude, as shown in the following equations. 
For the B Ki+i' decay. 



'2 kK 



d'B- 



m 



dsdu 



tb 



B 



277r3 



Re{WA'Al){X-u' 



+Re{WC'Al) - (A - u^) - 4:mj{2 + 2m 



+ |AiipmB[A - M + 2m^(2 + 2m\ - s) 



fBl 

dsdu 



fBl 

dsdu 



m 



Tb 



B 



2^3 



Re{WA'A*-){2meu) + Re{WC'Al){l - mj^){-2m 



-Re{WD'Al){-2mes) + |A^p(s - 2m^) | , 



m 



Tb 



B 



2^3 



Re{WA'Al){2miu) + Re{WC'Al){l - m^)(2m^) 



-Re{WD'Al){2mes) + \Ai,Ws -2m^g) | 



(52) 



(53) 



(54) 



with W = 
For the B 



dsdii 



2 kK* 



dsdu 



d^BT 



dsdu 



tb 



B 



decay, we have 



ReiWAI*) s{\ + u^) + AmjX 



-ReiWEV) [s(A + u^) - 4m^A] + |/|^ [s(A + u^) 

+Uu[Re{WAr) + Re{WBr) - Re{WEJ*) - ReiWFI*) + 2Re{IJ*) 
+ Re{WBr)(\ -u^ + 8ml, (s + 2m^) 
-Re{WFJ*)(^X -u^ + 8m%,{s - Amf 
+ \J\\\ -u^ + 8m|^.(s - m^)) 
-Re{WBK*){X - - m^. - s) 
+Re{WFK*)(^{X - u^){l-ml, - s) + Am^^X] 

-2Re{JK*){{X - u^){l - m]^, - s) + 2rh]X] 
X 



m 



K* 



Re{WCK*){X - n^) - Re{WGK*)(X - + 4m^(2 + 2m%, - s 



+ \K\'[X - r + 2m|(2 + 2m|^ 



ml. 



-X 



- Re{WHL*)s + \L\'s/2 + Re{WFL*) - Re{JL* 



-Re{WGL*){l - m\,) + Re{KL*){l - m\. 



(55) 



m 



Tb 



B 



mp 



2H^ 



m 



K* 



Im{WBni) ■ {X'^u{l - m|^. - s 



+Im{WCQl)X^u + Im{WFQl)X'^ 



-Im{WGni)X-^{l - ml,] 



+ |l]^|2(s-2m2) I, 



(56) 



Tb 



2^3 1 



m 



K* 



Im{WBn';) ■ (A~^n(l 



m 



K* 



^im{wcn';)x^u - im{WFn';)X'^ 



-Im{WGVt';)X-^{l-rhl,] 



'',?{s-2m]) 



(57) 



From the total double differential branching ratios, we can get the normalized forward- 
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backward asymmetries Afb PP 



+1 cfB{B~^K'-*h+e-) 



ds 



dsdcosO 



sign{cos6)dcos6 



dcosO 



(58) 



r+i 

J — 1 dsdcosd 

In the SM, the Afb vanishes in i? ^ Ki~^i~ decays as shown by Eq. flTTI) . since there is no term 
containing u with an odd power. The RPV effect via the squark exchange on AFBiB — >■ Ki~^i~) 
also vanishes for the same reason as shown by Eq.( l52i) . 

The total decay branching ratios of the pure leptonic Bg decays are calculated to be 

1 



B''\B, 
1 



1 + 



h 



SM\ 



+ 



+ 



\hsM\ 

1 



\h 



SM\ 



i?e(/i5A/A;)— + |Ap|' 



1 



2Re{hsMK:) + |Ad' 
1 1 



2m} 



rrii 



f 1 



1 



V2m| 



m 



Bs 



(59) 



From this equation, we can see that B{Bs ) could be enhanced very much by the 

s-channel RPV sneutrino exchange, but not by the t-channel squark exchange. 



3 Numerical results and analysis 

Now we are ready to present our numerical results and analysis. Firstly, we will show our 
estimations and compare them with the relevant experimental data. Then, we will consider 
the RPV effects to constrain the relevant RPV couplings from the recent experimental data. 
In addition, using the constrained parameter spaces, we will give the RPV SUSY predictions 
for B{Bs £~^i~) and Afb{B — > K^*'^i'^i~), which have not been well measured yet. 

For the form factors involving the B K^*^ transitions, we will use the recently light- 
cone QCD sum rules (LCSRs) results [23], which are renewed with radiative corrections to the 
leading twist wave functions and SU(3) breaking effects. For the dependence of the form 
factors, they can be parameterized in terms of simple formulae with two or three parameters. 
The form factors V, Aq and Ti are parameterized by 

For the form factors T3, /+ and /t, it is more appropriate to expand to second order around 
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the pole, yielding 

F{s) = + (61) 

1 — s/m^ (1 — s/mj^ 

where rh = rhfit for A2 and T3, and m = m/j for /_|_ and fx- The fit formula for Ai, T2 and /o is 

The form factor T3 can be obtained by T^{s) = ^"™a-* [T3(g) — T2(s)]. All the corresponding 
parameters for these form factors are collected in Table I. In the following numerical data 
analyses, the uncertainties induced by -F(0)[23] are also considered. 

Table I: Fit for form factors involving the B K^*^ transitions valid for general [23j . 





s) 


m 


^tot 


ri 


m| 


r2 




fit Eq. 




0.331 


0.041 


0.162 


5.412 


0.173 




m 




■*K 


0.358 


0.037 


0.161 


5.412 


0.198 




m 


/cf- 


■*K 


0.331 


0.041 






0.330 


37.46 


m 


yB- 


^K* 


0.411 


0.033 


0.923 


5.322 


-0.511 


49.40 


m 




0.374 


0.033 


1.364 


5.282 


-0.990 


36.78 


dnnD 


Af- 




0.292 


0.028 






0.290 


40.38 


m 


At 


*K* 


0.259 


0.027 


-0.084 




0.342 


52.00 


m 




*K* 


0.333 


0.028 


0.823 


5.322 


-0.491 


46.31 


ioD 


rpB- 


*K* 


0.333 


0.028 






0.333 


41.41 


m 


n^B- 
-'-3 




0.333 


0.028 


-0.036 




0.368 


48.10 


m 



The other input parameters and the experimental data are collected in Table II and III, 
respectively. In our numerical results, if not specified, we will study physics observables in 
the region s > O.lGeV"^, and use the input parameters and the experimental data which are 
varied randomly within la and 2cr variance, respectively. We assume that only one sfermion 
contributes at one time with a mass of 100 GeV. As for other values of the sfermion masses, 
the bounds on the couplings in this paper can be easily obtained by scaling them with factor 

f2 = ( \2 
J — viOOGeV/ • 
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Table II: Default values of the input parameters and the itlo" error bars for the sensitive param- 
eters used in our numerical calculations. 

me, = 5.370 GeV, m^^ = 5.279 GeV, = 5.279 GeV, mw = 80.425 GeV, 

mK± = 0.494 GeV, trko = 0.498 GeV, mK'± = 0.892 GeV, uik-o = 0.896 GeV, 



fribimb) = (4.20 ± 0.07) GeV, ms{2GeV) = (0.095 ± 0.025) GeV, 




mu{2GeV) = 0.0015 ~ 0.003 GeV, md{2GeV) = 0.003 ~ 0.007 GeV, 




me = 0.511 X 10-3 GeV, = 0.106 GeV, mt,poie = 174.2 ± 3.3 GeV 


m 


tb, = (1.466 ±0.059) ps, r^, = (1.530 ± 0.009) ps, r^^ = (1.638 ± 0.011) ps. PU 


\Vtb\ ^ 0.99910, \Vts\ = 0.0416lt°:S. 




sin2^v^' = 0.22306, = 1/137. 


[27\ 


fn, = 0.230 ± 0.030 GeV [2H] 



Table III: The SM predictions and the experimental data for B 
K^*k+£- and Bs t+r [D El H EH [29] . 





SM prediction value 


Experimental data 


B{B - 




0.610x10-6 


(0.56ll°:°66) X 10-6 


B{B - 


^ Ke+e-) 


0.610 X 10-6 


(0.380l°r7) X 10-6 


B{B - 




1.27 X 10-6 


(1.44 ±0.23) X 10-6 


B{B - 


^ K*e+e~) 


1.29 X 10-6 


(1.25 ±0.27) X 10-6 


B{Bs- 




3.72 X 10-9 


< 1.0 X 10-^ (95% C.L.) 


B{Bs- 


e+e-) 


8.70 X 10-1^ 


< 5.4 X 10-5 (90% C.L.) 



The branching ratios in the SM estimated with the central values of the input parameters 
are presented in Table III, and the relevant experimental data P, El HI E3 EHj are listed for 
comparison. From Table III, we can find that the branching ratios of the semileptonic decays 
roughly consistent with the SM predictions. So, there are still windows, however limited, for 
RPV contributions. 

We now turn to the RPV effects. There are six RPV coupling products contributing to 
four B — i> K^*k~^£~ and two Bg i^£~ decay modes. We use the experimental data of the 
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branching ratios listed in Table III to constrain the relevant RPV coupling products. 

For the RPV couplings \iii\%2 A*^]^A'j23, since their RPV weak phases are found to have 
very small contributions to the physical observables, we take their RPV weak phases to be free, 
and only give the upper limits for their modulus which are listed in Table IV. In FiglHl we 
present our bounds on the other four RPV coupling products. 




Figure 3: The allowed parameter spaces for the relevant RPV coupling products constrained 
by the measurements listed in Table IV, and the RPV weak phase is given in degree. 

From Figjni we find that every RPV weak phase are not constrained so much, but the 
modulus of the relevant RPV coupling products are tightly upper limited. The upper limits for 
the relevant RPV coupling products by -B ^ K^*H~^i~ and Bg decays are summarized 

in Table IV. For comparison, the existing bounds on these quadric coupling products [221 EO] 
are also listed. From Table IV, one can find that our bounds are more restrict than the previous 
ones. It's also noted that the previous bounds of |AiiiA'j32| and |A*^iA^23l obtained at the 
Mgut scale in the RPV mSUGRA model [30]. 

Using the constrained parameter spaces shown in Table IV and FigOl we can predict the 
RPV effects on the other quantities which have not been well measured yet in these processes. 
We perform a scan over the input parameters and the new constrained RPV coupling spaces 
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Table IV: Bounds for the relevant RPV coupling products by B — > K^*H~^i~ and 
Bg decays for 100 GeV sfermions, and previous bounds are listed for com- 

parison. 



Coupling 


S 


Bounds [Processes] 


Previous bounds [Processes] 


1 \' \'* 

I^ill^i23 
I^2j3'^2i2 
Ai22A-32 

I^i22^i23 




^ /I 7 v 1 n— 5 [R ^ 1 

^ 4. ( X iU [JD — >■ l\ ^ ' e 6 J 

< 2.3 X 10-5 [B K^*^e+e-] 

< 2.3 X 10-5 [B K^*h+e~] 

< 4.6 X 10-5 

< 1.8 X 10-5 

< 1.7x10-5 ^Xl-^ 


^ y.D X iu [JD — > ixe 6 J [zz] 

< 1.5 X 10-3 [30] 

< 1.9 X 10-^ [30] 

<3.9xl0~3 [g^^^+^-j 
<9.7xl0-5 [B^Kfi+fj.-] 1^ 

< 2.7 X 10-5 [B JsT/iV] [22] 

< 2.7 X 10-5 [B Kji+ii-] [22] 



to get the allowed ranges for B{Bs and Afb{B K'-*H+i~) with the different RPV 

coupling products. Our numerical results are summarized in Table V. 



Table V: The theoretical predictions for B{Bs and Afb{B ^ K^*^i+(.-) in the SM and the 

RPV SUSY. The RPV SUSY predictions are obtained by the constrained regions of the different 
RPV coupling products. The index g = 1 and 2 for i = e and ;U, respectively. 





SM value 




-^199^^2 


\gg\23 


B{B,^e+e~) 


[5.9,11.2] X lO-i'' 


[5.3, 17.0] X 10-1^ 


< 2.2 X lO-'^ 


< 2.2 X 10-^ 


Afb{B ^ Ke+e-) 








[-4.1,5.0] X 10-5 


[-4.0,5.0] X 10-5 


Afb{B ^ K*e+e-) 


[0.12,0.18] 


[0.00,0.29] 


[0.12,0.18] 


[0.13,0.18] 


B{Bs^fi+fi') 


[2.7,4.8] X 10-9 


[2.1,6.8] X 10-9 


< 1.0 X lO-'^ 


< 1.0 X lO-'^ 


Afb{B ^Kfi+fi-) 








[-7.0,6.9] X 10-3 


[-6.8,7.0] X 10-3 


Afb{B ^ K*i^+^r) 


[0.13,0.18] 


[0.02,0.76] 


[0.13,0.18] 


[0.13,0.18] 



From the Table V, we can find some salient features of the numerical results. 

• I. As shown by Figj2](b), the contributions of A'^j3A'i*2 and A2j3A2*2 to B{Bs — > e^e~) and 
B{Bs fi^fi~), respectively, arise from t-channel squark exchange . After Fierz transfor- 
mation, the effective Hamiltonian due to the t-channel squark exchange is proportional 
to (^s'y'^'Pjibj {I'^^Pitj, which contribution to B{Bs is suppressed by mj/rn^ due 
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to helicity suppression. Therefore B{Bs e~^e~ , will not be enhanced so much by 

the t-channel squark exchanging RPV contributions. However, the effective Hamiltonian 
of s-channel sneutrino exchange would be ± 75)6] (^(1 =F 75)^)) which contributions 



are not suppressed by m\/m\. So that, both B{Bs — > e+e ) and E{Ba — > ) could 
be enhanced to order 10^'' by AjiiA^32(A,*^]^A'j23) and Ai22A^32(A*22A^23)5 respectively. 

• II. Both Afb{B Ke+e-) and ApsiB Kn+n') are zero in the SM. The RPV 
contributions to the asymmetries due to squark exchange are also zero, while the sneutino 
exchange RPV contributions are too small to be accessible at LHC 

• III. It is interesting to note that the RPV squark exchange contributions have signifi- 
cant impacts on Afb{B K*e~^e~) and Afb{B K*fi~^fi~). However, the sneutrino 
exchange have negligible effects on Afb{B K*e'^e~) and Afb{B — > K*fi~^fi~). 

Recently, the BABAR Collaboration [1] has measured 

Afb{B+ ^ K+tn{s>^.iGeV^) = 0.151°:^^ ± 0.08, (63) 
Afb{B ^ K*tn{s>o.iGeV-) > 0.55 (95%aL.), (64) 

and Belle Collaboration has measured the integrated forward-backward asymmetries Afb [S] 

Afb{B^ K+tr) = 0.10 ± 0.14 ± 0.01, (65) 
AFBiB K*tr ) = 0.50 ± 0.15 ± 0.02, (66) 

where ^j?^ is slightly different from Afb, and defined by 

AFBiB ii)- . (67) 

J dcostids 

We find that the SM predictions for Afb{B — > K*i~^i~) are consistent with both the measure- 
ments within error bars. 

In FigsJUini we present correlations between the physical observable B, Afb and the param- 
eter spaces of the different RPV coupling products by the three-dimensional scatter plots and 
the two-dimensional scatter plots, respectively. The dilepton invariant mass distribution and 
the normalized forward-backward asymmetry are given with VMD contribution excluded in 
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terms of dB/ds and dAps/ds, and included in dB' /ds and dA'p^/ds, respectively. From FigsHl- 
[Hl one can find the correlations of these physical observables with RPV coupling products. In 
FigslMSl since the influences of the RPV weak phases are very small, we take the RPV weak 
phases randomly varied in [—it, it] and only give the change trends of the physical observables 
with the relative modulus in the two-dimensional scatter plots. 

At first, we will discuss plots of FigJUin detail. The three-dimensional scatter plot Fig|l](a) 
shows Afb{B K*e~^e~) correlated with |A'i,j3A']*2l and its phase (p^^ . We also give projections 
on three vertical planes, where the |A'ij3A']*2|-0^p plane displays the constrained regions of 
Xii-s\ii2 as the first plot of FigOl It's shown that Afb{B K*e^e~) is decreasing with |A'ij3A']*2l 
on the Afb{B K*e~^e~)-\\'n-^\ii2\ plane. From the Afb{B — > K*e~^e~)-(j)^j^ plane, we can 
see that Afb{B — > K*e~^e~) is decreasing with 10^^ |. The recent measurement of Afb{B — >• 
K*i~^i~) favors \^'ii3^ii2\^ (2 ~ 4) x 10"^ with small (f)^^. However Bs e^e~ is remained 
inaccessible at LHC From plots Fig|l](c-h), we can find that the theoretical uncertainties in the 
calculation of the SM contributions are still very large, nevertheless, for dB{B — > K^*^ e~) / ds 
and dApB^B — > K*e~^e~)/ds, the A'^j3A'^*2 contributions are distinguishable from the hadronic 
uncertainties. 

Fig|5]and Fig|6]are the contributions of AjiiA^32 and A*]^^A^23) respectively. It is interesting to 
note that the t-channel RPV sneutrino exchange could enhance B{Bs e^e^) by about 6 orders 
of magnitude (^Figj^b) and Figl6](b)^. Moreover, they could give distinguishable contributions 
to dB{B Ke~^e~)/ds in the high s region. For other observables, their contributions are 
indistinguishable from hadronic uncertainties. 

The RPV squark exchange contributions to i? ^ K^*^ii'^^~ and Bg 1^^ are presented 
in FigJTl Such contributions could give large Afb{B K*fi^fi~) in favor of Belle and BABAR 
measurements, however, small corrections to B{Bs 

The RPV sneutrino exchange contributions to 5 — K*fi~^fi~ and Bg fi~^fi~ are displayed 
in FiglH] and FiglHl respectively. From plots in these two figures, we find B{Bs —>■ fJ'~^fi~) 
is very sensitive to such contributions. Within the constrained parameter space for Aj22A^32 
and A*22A'j23, RPV contributions can enhance B{Bs f^^fJ' ) by two orders, which could be 
accessible at LHC and Tevatron in the forthcoming years. We also note that these contributions 
to other observables are small, regarding to the large theoretical uncertainties. 
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It is worth to note that, in the BABAR measurement [T] , the value of Afb{B K*£~^£~) 
with 95% C.L. lower limit is slightly above the SM prediction in low s region. We have found 
that the RPV couplings A'^^gA'^^jg and A2j3A2*2 could enhance Afb{B K*£~^i~) to accommodate 
the possible discrepancy, which are shown by the following Fig. 4(g) and Fig. 7(g), respectively. 

Generally, the predictions of B decays suffer from many theoretical uncertainties. Recently, 
beyond naive factorization, B K* have been investigated with the QCD factorization [31] 
framework in [321 and then studied with Soft Collinear Effective Theory (SCET) [M] in 
Ref . [35] . However, the power corrections are found to be associated with large theoretical 
uncertainties, and we are conservative to the naive factorization in this paper. 

To probe new physics effects, it would be very useful to measure correlative observables, 
for example, dApB^B — > K* fi~^ fi~)/ds, dB{B K* fi~^ fi~) / ds and B{Bs — > /i^/i~), since 
correlations among these observables could provide very strict bound on new physics models. 
At present, one may have to wait for the error bars in the measurements of these observables to 
come down and more channels to be measured. With the operation of B factory experiments, 
large amounts of experimental data on hadronic B meson decays are being collected, and 
measurements of previously known observable will become more precise. From the comparison 
of our predictions in FigsHHO] with the near future experiments, one will obtain more stringent 
bounds on the products of RPV couplings. On the other hand, the RPV SUSY predictions of 
other decays will become more precise by the more stringent bounds on the RPV couplings. 

4 Conclusions 

In conclusion, the decays B K^*^£'^i~ and Bs — * are very promising means to probe 
effects of new physics scenarios. In this paper, we have studied the B K^*k+£- and B, 

decays in the RPV SUSY model. We have obtained fairly constrained parameter spaces 
of the RPV coupling products from the present experimental data of B ^ K^*H'^i~ and Bg — 
£~^i~ decays, and found these constraints are stronger than the existing ones, which may be 
useful for further studies of the RPV SUSY phenomenology. Furthermore, using the constrained 
parameter spaces, we have presented the RPV effects on the forward-backward asymmetries 
of -B — > K^*H~^i~ and the branching ratios of the pure leptonic Bg decays. Our results of 
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Afb{B K^*^i~^£~) agree with the recent experimental data. It is shown that B{Bs 
could be enhanced several orders by the RPV couplings from the sneutrino exchange. Since we 
have poorly experimental information about the pure leptonic decay Bs —* , further refined 
measurements of B{Bs can further restrict the constrained space of the four RPV 

couplings from the sneutrino exchange. We have also compared the SM predictions with the 
RPV predictions about dilepton invariant mass spectra and the normalized forward-backward 
asymmetries in i? ^ K^*H~^i~ decays. We have found that dB{B — > Ke~^e~)/ ds could be mildly 
decreased by the RPV coupling A'^j3A']*2, and the A2j3A2*2 contributions to dB{B Kji'^ n~) / ds 
are indistinguishable from the hadronic uncertainties. The other four RPV couplings due to 
the sneutrino exchange have distinguishable effects on dB{B K£'^£~)/ds at high s. For the 
B — ^ K*i'^i~ decays, the squark exchange contributions have distinguishable effects on dB/ds 
and dApB/ds, but the contributions from the sneutrino exchange have no apparent effects on 
them. The results in this paper could be useful for probing RPV SUSY effects and will correlate 
strongly with searches for direct RPV signals at LHC in the forthcoming year. 
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Figure 4: The effects of RPV coupling A'^j3A']*2 due to the squark exchange in i? ^ K^*^e~^e~ 
and Bs — e~^e~ decays. The primed observables are given with ip{nS) VMD contributions, and 
-07 denote s. 
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Figure 5: The effects of RPV coupling AjiiA^32 due to the sneutrino exchange in i? — > K^*^e~^e~ 
and Bs — * e~^e~ decays. 
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Figure 6: The effects of RPV coupling \*ii\'i23 due to the sneutrino exchange in i? — > K^*^e~^e~ 
and Bs — * e~^e~ decays. 
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Figure 7: The effects of RPV coupling \'2i'i^'2i2 due to the squark exchange in S — > K^*^n~^n 
and Bg decays. 
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(9) 



Figure 9: The effects of RPV coupling \I22K2?, due to the sneutrino exchange in S — > K^*^n~^n 
and Bg decays. 
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